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The effect of finite compressibility of the electron gas is taken into account in the calculation of the propa- 
gation constant and the power carried by a VLF wave propagating in the ionosphere. 

The principal result is that the electromagnetic wave ("whistler" mode) and a quasi-acoustic or com- 
pressional wave are coupled near the critical, or maximum, angle of the wave normal with respect to the direc- 
tion of the geomagnetic field. Near the critical angle, the direction of power flow is no longer const rained to 
lie within a narrow cone of angles about the geomagnetic field, but energy can propagate in a quasi-acoustic 
mode at all angles up to the critical angle of the wave normal. Because of collisional (lamping, the effects of 
finite compressibility do not appear to be detectable unless the electron temperature in the outer ionosphere 
exceeds 10 4 °K. 

1. Introduction 

Recent experiments in communication with satellites and space probes have stimulated 
considerable interest in studies of propagation of electromagnetic waxes in ionized gases. This 
interest has not been confined to the higher radiofrequencies, but extends to the VLF band as 
well. Although much of the basic phenomenon of VLF propagation within the ionosphere 
has been satisfactorily explained by Storey [1953] and by more recent workers (for example, 
[Helliwell and Morgan, 1959; Smith, Helliwell, and Yabroff, 1960; Yabroff, 1961]), in connec- 
tion with the study of atmospheric "whistlers," the possibility of performing satellite experi- 
ments [Leiphart et al., L962; Cain et ah, 1961] at VLF, utilizing antennas located within the 
ionosphere, affords the* opportunity of measuring other effects which may not be measured 
directly from the ground. 

One of the effects which lias not previously been accounted for in the theory of ionospheric 
propagation is the finite compressibility of the electron pis. The effects of finite compressi- 
bility may not be apparent at frequencies above the gyro frequency, because of the great dis- 
crepancy in the propagation velocities of electromagnetic and acoustic waves. At VLF, 
however, the great phase retardation of the electromagnetic wave which can occur appears to 
provide favorable conditions for interaction of the electromagnetic wave with an acoustic, or 
compressional wave. It is the purpose of this paper to invest igate the effects of finite compressi- 
bility on the propagation of small-signal, plane waves at VLF frequencies in a uniform 
ionosphere. 

2. Magneto-ionic Theory for an Electron Gas 

In considering propagation at frequencies of several kilohertz and higher, the motions of 
heavy ions may be ignored, since the ion plasma frequency and ionic gyrofrequency in the 
earth's atmosphere are less than a few hundred hertz. 

The field equations for a time-harmonic (e jwt ) electromagnetic field in an electron gas can 
then be written [Oster, 1960] in linearized form for fields of small amplitude (fx II r f<^<CB ) y 
and for zero average drift velocity of the electron gas, 

VXE=-jco Mo H (1) 

VXH=jue E-Ne\ (2) 



1 Pan of the results presented in this paper were given in a paper presented at the IEEE Symposium on Antennas and Propagation in Boul- 
der, Colo., July 1963. 
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V-v=-^ (3) 

Vp=-(JG)+v)mNv-Ne(E+\XB ) (4) 

where 

E, H represent the time-varying electromagnetic field intensities, 
B is the flux density of the geomagnetic field, 
N is the average electron density, 

v is the electron velocity, averaged over a small element of volume whose radius is the Debye 
length, 

e is electronic charge, 
m is electron mass, 

v is the electron collision frequency, 
P is the time-averaged pressure of the electron gas, 

p is the time-varying component of pressure, 

7 is the ratio of specific heats, equal to 5/3 for a gas of free electrons. 

Equations (1) and (2) are just Maxwell's equations, taking into account convection 
currents due to electron motion; (3) is the combined equation of continuity and state, and 
(4) is the momentum equation. 

Equations (3) and (4) can be combined to give 

-^WV=-(jco+^)miVv-^(E+vXBo) (5) 

which now gives the behavior of the electron plasma, expressed in terms of the velocity, in 
response to the electromagnetic forces acting upon it. 

Now to solve the field equations (1), (2), and (5), consider propagation of a plane wave. 
Expressing the field vectors in terms of a Cartesian coordinate system, and taking the £-axis 
along the direction of propagation, it is seen that the spatial dependence of all time-varying 
quantities is given by e~ 3 ^ z , where /3 is the propagation constant, which may be complex. 

The two Maxwell equations, (1) and (2), may be combined to give a wave equation, 

V X V X E-kSE=ja»oNeY (6) 

where & =Wmo€o is the wave number in free space. In view of the assumed plane-wave spatial 
dependence of the fields, (6) becomes 

(^-^ 2 -w)E=-ico Mo iV6V. (7) 

It is convenient to write (7) in matrix form, observing that (V-E)^^ 2 ^, obtaining 

"*8— j8 2 0" 

kl-p 2 

o o &|L 



. e 
J com 



•E=«v (8) 



where k„= — — = coa/* eo is the square of the wave number corresponding to the electron 
plasma frequency. Similarly, by taking the coordinate system such that B lies in the x-z 
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plane, (5) can be written as 



E= 



. com 



1-3- 



3% 

CO 



*>L 



i-i 



. co T 



-J- 



CO 



2 

CO" 



(9) 



.OOm -r 

J e 



where 



_rA : 



u2= ~aj~~ ls the square of the acoustic velocity in the medium, 

co L =co cos 0, 
co T =co G sin 0, 

6 

coo= — B 0y the electron radian gyrofrequency, 

m ° ^ J 

6 is the angle between the direction of the wave normal and the geomagnetic field. 



ykT 

m ' 



wnere 



The acoustic velocity is directly related to the temperature of the electron gas, u 2 -- 

Tis in degrees Kelvin, and k is the Boltzmann constant. 

Combining (8) and (9) and eliminating E gives a determinantal equation for the propa- 
gation constant, 



det 



(l-J-j)(«-«-« 







0) 







CO 



=0. 



(10) 



CO 

i*« (i-3 v -)%-^-K 

CO \ CO/ CO 

It is the solution of (10) for values of (3 which is now of concern. 

3. Perturbation Solution of the Determinantal Equation 

The determinantal (10) yields a cubic equation in 0\ with complex coefficients. An analytic 
solution of the equation is difficult except for some special cases, when the direction of propa- 
gation is either along, or transverse to the direction of the geomagnetic held. For the trans- 
verse case, there are no solutions which yield Re (/3 2 )>0, or a propagating wave. In the longi- 
tudinal case, two such solutions occur. One of these values of /5 2 corresponds to whistler mode 
propagation, and the other to a longitudinal acoustic wave. 

A perturbation solution to the determinantal equation may be obtained using the solutions 
for longitudinal propagation as the zero-order solution. The first perturbation of the equation 
will be shown to give a propagation constant for the electromagnetic wave which is recognizable 
as the quasi-longitudinal approximation of magneto-ionic theory for an incompressible plasma, 
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plus another propagation constant for an acoustic wave. Application of the perturbation 
method to obtain a second order correction is not fruitful, since the perturbation parameter is 
not small for all angles 0. 

It is useful to rewrite (10) as an operator equation 



-(l-jZ)(l-\)-X -jY L (l-\) 

Sx-v= jY L (l-\) (i-jZ)(i-\)-X 

jY T 

where, in the usual notation for the magneto-ionic theory, 



X- 



"n l^n 



0>L 





-jF r (l-X) 

l-jZ-V\-Xj 



v=0 



(11) 



Y T = 



0) T 



and, in addition, 



V- 



w 



and X is the square of the refractive index. It can be seen that X is an eigenvalue of an operator 
S by separating the terms containing X from S\, 

(12) 



(13) 





Sx-v= 


=[S -\ix\-y-- 


=0. 


The operators S c and n are, in matrix form, 






-p-x 


-JY L 


" 


o = 


JY L 


P-X 


-JY T 







jY T 


p-x. 




- p 


-JY L 


' 


M= 


JY L 


P 


-JY T 










V 



(14) 



where P=l—jZ. 

The eigenvalue equation (12) can be put into a more conventional form by a coordinate 
transformation diagonalizing fi. This can be accomplished by defining 



which gives the eigenvalue equation 

where 

1 



>o * M 



PX 



P 2 -Y 2 L 

XY±_ 
J P 2_yi 

Y L Y T 



x=/*.v, 

r-x-xx-o 

J JJ2 y-2 

_px_ 

P 2 -Yl 
. PY T 

J 7)2 TZ2 
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(15) 
(16) 



XY L Y T 

V(P 2 -Yl) 

. PXY T 



F(P 2 -F|) 



P-X 



PY% 



V V(P 2 -Yl). 



(17) 



The operator T is now written as the sum of a zero-order term and a perturbation term, 



T=T Q +T, 



(18) 



where 



PX 



XY L 



■ XY L 
J >-F£ 

PX 







P-X 

V 



(19) 



and 



2i= 



r, 



jp 



xy, ■ 

V 
.PX 

PY T 
V 



(20) 



Thus To corresponds to the quasi-longitudinal a pproximati on of magneto-ionic theory, and the 
expansion parameter of the perturbation is Yr/^lYl—P 2 . 

The eigenvalues and corresponding normalized eigenvectors of 7',, are 



\!°'=i- 



X 



P-Y 




W2 V2 / 

■ VV2 

xf = (0, 



* 



o, 



1). 



(21a) 
(21b) 
(21c) 



In calculating the eigenvectors, the damping has been neglected, by taking P=l. This is 
justifiable for the case of small losses since the eigenvectors are to be used only in power cal- 
culations, in which a vector is multiplied by its complex conjugate in order to find real-valued 
quantities. In such calculations, the error incurred in neglecting Z in the eigenvectors will be 
of order Z 2 in the power calculations, and can therefore be neglected for Z 2 <<Cl. 

The first-order approximation to the eigenvalues is now calculated from the zero-order 
eigenvectors and the operator T, 



r-x? 



These are 



xsf ) 



At " 


X 


(Q) . v(0) 


x?» 


=1- 


X 
P-Y L 


*?>= 


=1- 


X 

P+Y L 


P- 


-X 


PY% 



V V(P 2 -Yl) 



(22) 

(23a) 
(23b) 
(23c) 



Since V is ordinarily a very small number, only terms containing negative powers of V have 
been retained in the expression for X£ l) . 
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For Z<<0, the complex eigenvalues can be expanded to first power in Z, giving 

X ' )=1 -T^( 1+ ^> ^ 

X^l-^l+j^} (24b) 



1-X-jZ Y\ 

V V{\-Yl) 



( 1+ fT=r} (24c) 



It can be seen that the first two eigenvalues still agree with the results of the QL approxi- 
mation of magneto-ionic theory, in which only one of these eigenvalues has a positive real 
part for X>F L >1. The third eigenvalue has a positive real part for F|/(F£-— 1)>X—1, 
which indicates that a quasi-acoustic mode can propagate over a narrow cone of wave-normal 
angles near the critical angle for whistler mode propagation, which is given approximately by 

c ^arccos y 

The damping of either the electromagnetic whistler mode or of the quasi-acoustic mode be- 
comes large near the critical angle, thus the acoustic mode may not be observed due to exces- 
sive damping if the cone of the angles along which it may propagate is too small. The width 
of this cone can be shown to be approximately given by 



cos20 mla — cos20 max ~ 



2 

X 



Thus the acoustic wave with low damping can probably be observed only in the outer iono- 
sphere, above the height of maximum electron density in the F2 layer. 

It is now of interest to examine the first-order approximation to the eigenvectors. These 
are found from perturbation theory by writing the eigenvalue equation as 

[T-A].x=[T -X (0) ]-x (0) +S[^ (25) 

where 5 is a dummy perturbation parameter. Since [ T — X (0) ] .x (0) = 0, x (1) can be found by 
setting 

[r -X (0) ]-(x (1) -x w ) = -[7 T 1 -(X (1) -X (0) )]-x (0) (26) 

and solving for the vectors x (1) subject to the restriction 

(x (1) -x (0) ).x (0) =0. (27) 

This restriction enables x (1) to be specified uniquely. The new set of vectors which are found 
by this procedure are, again neglecting Z in solving for the eigenvectors, 

*r=( l -i (i-VHi-r?-X)+AY} <28,) 

*f'-(u- (i-r)( 1+ £-X)-^ ) (»> 

(1) / XY t • XY t -,\ zoo \ 

* 3 -\~(l-X)(l-Yl)' J {l-X)(l-Yiy l f [Zbc) 

4. Poynting Vector 

In order to find the ray direction, or direction of energy transport, it is necessary to evaluate 
the power flow from the generalized Poynting vector. Thus for a wave having both electro- 
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magnetic and acoustic energy, the power density is given by 

P=|Re[EXH+pv acougtlc ]. (29) 

For plane waves with the wave normal along the z coordinate, this expression can be written as 

P=f Re /^[E.Ea g -^+(^yZF^a 2 iy (30) 

To find the quantities of interest, v and E from x, it is necessary to invert the transforma- 
tion (15), finding 

v=m _1 -x. (31) 

This form is appropriate to find v z for use in the power calculation (30), but is not appropriate 
to find E. The electromagnetic energy term should not depend on the longitudinal velocity 
component, so we find that part of the electric field which contributes to the electromagnetic 
power flow by employing a projection operator, 



Q- 



Then, for use in the calculation of (30), 



1 





0" 





1 











()_ 



^ .com - r , ~ 



(32) 



(33) 



where L is the operator defined by (9). 
One finds then that 

W=-j^ (l,-j,Y T /(l-Y L )), 
vV=Y T /(X(V-l+Y L )), 
and, for Z«l so that X/(X— 1)~1, 

^=jfj^Yi (W, Y T /(1-Y L )), 

From these vectors, the power expressions are easily calculated, 
~^\ e ) LV + X(V-l+Y L y 



(34a) 



£0jU 






(34b) 



(35a) 



(35b) 



The ray directions, or directions of power flow, deviate from the wave-normal directions by 
angles 

&=arctan (pr\ i=l, 3. (36) 

Because the ratio X/Vis a very large number, it is seen that the ray direction of the quasi- 
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acoustic wave, P 3 , is nearly along the wave normal, whereas that of the quasi-electromagnetic 
wave, Pi, becomes very nearly in the direction of the geomagnetic field for angles near the 
critical angle, if F>>1. 

The wave impedance of the quasi-electromagnetic and the acoustic waves can also be 
determined from the power calculation. Defining a wave impedance in the direction of the 
wave normal by 

the impedance is found to be 



Z 2 = 



WMo 



l+^J — =- 

E-E-E Z E 2 



(37) 



In the quasi-electromagnetic wave, the ratio of acoustic to electromagnetic energy is 
small, and the impedance becomes approximately 



z lt 



"ft ' 



(38) 



For the quasi-acoustic wave, the acoustic energy is many times greater than the electromag- 
netic energy, and the wave impedance is correspondingly greater, 



Zr 



co/z /com 



XV — — 



v 2 v z 



Es'Ea — EszEa: 



(39) 



As a result of the disparity in wave impedances, an antenna which is designed to radiate an 
electromagnetic wave efficiently will not couple greatly to the quasi-acoustic wave, and in order 
to observe acoustic wave propagation in the ionosphere, it may be necessary to employ special 
antennas. 

5. Some Numerical Results 

In figures 1 and 2 the refractive index is given for a lossless ionosphere (Z=0) as a function 
of wave-normal angle for parameters X=50, F=7.8, F=0.002 2 , as obtained by a numerical 
solution of the exact determinantal (10). The refractive indices obtained by means of the 
first-order perturbation solution are also given in figure 2 for comparison. The longitudinal 
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Figure 1. Section of refractive index surface for X = 50, Y=7.8, Y=0.002. 

The angle 8 is the inclination of the wave normal with respect to the geomagnetic field. The arrows along the curve indicate directions of 
power flow for waves having various wave-normal angles 0. 



2 These parameter values correspond roughly to a model of the earth's ionosphere at a distance of slightly more than one earth radius above the 
earth's surface. The value of V corresponds to an electron temperature of about 7X10" °K. These values were not chosen to give a realistic 
approximation to any actual conditions so much as they were chosen to give results which could be illustrated graphically. A more realistic 
set of values would have X~400— 7000 and V~5X(10~ 5 — 10 -7 ), yielding a larger critical angle and a narrower cone of angles in which the quasi- 
acoustic wave can propagate. 
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Figure 2. Refractive index as a function of tan 6, 
indicating the transition from a quasi-electromagnetic 
to a quasi-acoustic wave. 

The first-order perturbation solution is compared with the exact 
solution of the determinental equation for X=50, Y=7.8, F=0.002. 



acoustic wave is not shown in these figures. It is confined to wave-normal angles very near 
the direction of the geomagnetic field and does not couple with the electromagnetic wave. 

The curve of figure 1 shows the refractive index and the direction of power flow, or ray 
direction, versus wave-normal angle with respect to the direction of the geomagnetic field. 
Up to the angle of transition from the quasi-electromagnetic to the quasi-acoustic wave, the 
refractive index is similar to that for the whistler mode in a cold or incompressible ionosphere, 
and the ray direction is approximately normal to the plane tangent to the refractive index- 
surface, as for the incompressible case. For larger angles, the wave becomes primarily acoustic 
in nature; the refractive index increases very rapidly with the wave normal angle, and the 
ray direction is nearly along the wave normal. 

From the approximate form (24c) for the square of refractive index, it is seen that the 
quasi-acoustic wave exists from an angle d A to C , where the angles are related by 



(40) 



6. Damping Effects 



Damping effects can be accounted for by collisions between electrons and heavy particles, 
principally positive ions and neutral atoms. The scattering cross section presented by charged 
ions to electrons is relatively very large at low electron temperatures, due to the long range of the 
Coulomb force. However, as the temperature, or mean kinetic energy of electrons increases, 
the scattering cross section of ions decreases, until at approximately 1.8 X10 5 °K the scattering 
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cross section presented by the Coulomb potential is no greater than is that of short range forces. 
The result is that the collision frequency in the ionosphere can be extremely low if the electron 
temperature is of the order of 10 5 °K. For instance, if X=50, T=10 5 °K, and it is assumed 
that the density of positive ions is equal to that of electrons, while the density of all heavy 
particles is 10 times that of electrons, then the collision frequency, estimated by the formulas 
given by Cowling [1945] would be approximately 8 X 10~ 4 , giving Z at 18 kHz equal to 4.4 X 10~ 8 . 
The damping of the quasi-acoustic wave will be greater than that of the quasi-electromag- 
netic wave, particularly near the region of coupling of the two waves. The ratio of attenuation 
constant to propagation constant for each of the waves is approximately 

tan * 1= Re^F^2F^l ? (41a) 

tan * 3 ~ Re ^ ^2 (Y L -l^l-(l-X)(l-Yl)IY* T (41b) 

Near the transition angle, the denominator of the term on the right of (41b) is small, giving 
a larger damping term than that of the quasi-electromagnetic wave. At the transition angle, 
damping will be greater, the smaller the value of V. For very small values of V, corresponding 
to electron temperatures in the range 10 3 °K or less, the damping will probably be so great 
as to prevent observation of a quasi-acoustic wave at wave-normal angles near the transition 
angle in any part of the ionosphere. 

7. Conclusions 

From the analysis presented here it appears that quasi-acoustic waves can be propagated 
within the ionosphere at VLF frequencies provided the electron temperature is sufficiently 
high (10 4 -10 5 °K). These waves will be coupled with the propagating electromagnetic waves, 
the so-called "whistler" mode of magneto-ionic theory, near the critical angle between the 
wave normal and the direction of the geomagnetic field. 

The effect of losses due to collisions between electrons and heavy particles is important, 
and in general the quasi-acoustic wave is subject to greater damping than is the electromag- 
netic wave. The damping of the quasi-acoustic wave increases with increasing electron density 
and decreases with increasing electron temperature. Damping prevents observation of the 
quasi-acoustic mode in the lower part of the ionosphere. 

It also appears unlikely that waves propagating in a natural whistler mode which couples 
in the outer ionosphere to an acoustic wave will be able to be detected at the earth's surface, 
although this conjecture has not been confirmed by ray- tracing calculations. However, 
only waves which have wave normal directions near the critical angle with respect to the 
geomagnetic field will be appreciably coupled, and such waves may either be "lost," leaving 
the earth's field, or may be highly attenuated in the lower region of the ionosphere. 

Acoustic waves at VLF may also be directly excited by radio transmitters located on earth 
satellites or by large-scale disturbances, such as nuclear explosions in the upper atmosphere. 
The possibility of excitation of such waves should be taken into account in the design of VLF 
satellite transmitting antennas, as the effect on the antenna impedance could be appreciable 
for certain orientations with respect to the geomagnetic field. 
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